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Abstract 

We present a study of the heavy-light spectrum and of the D- and £>-meson decay 
constants. The results were obtained in the quenched approximation, by using the non- 
perturbatively improved Clover lattice action at (3 = 6.2, with a sample of 100 configurations, 
on a 24 3 x 64 lattice. After a careful analysis of the systematic errors present in the extrac- 
tion of the physical results, by assuming quite conservative discretization errors, we find 
f Da = 231 ± 121? MeV, f D = 211 ± 14l° 2 MeV, f D Jf D = 1.10(2), f Bs = 204 ± 16+f MeV, 
fs = 179 db 18j]g 6 MeV, /b s //b = 1-14(3)15- ^ ur results, which have smaller discretization 
errors than many previous estimates at fixed value of the lattice spacing a, support a large 
value of fs in the quenched approximation. 



PACS numbers: 12.38.Gc,12.39.Hg,14.40.NdLb,13.25.Hw,13.30.Ce. 
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1 Introduction 



In this paper, we present the results of a lattice calculation of physical quantities of phe- 
nomenological interest for heavy quarks, such as their mass spectrum and decay constants. 

In order to reduce the systematic errors, we have performed calculations using the most 
recent developments in the lattice approach, namely: 

1. The non-perturbatively improved lattice Clover action [J3J, which we denote as "Alpha 
action" |2[ ||, || (see also ||), with the coefficient of the chromomagnetic operator 
computed in Ref. Q; 

2. Non-perturbatively improved vector and axial-vector currents, the renormalization co- 
efficients of which have been computed, using the Ward Identities method ^, [||, 0], in 
Refe. [|||; 

The use of non-perturbatively improved actions and operators allows us to reduce the dis- 
cretization errors to 0(a 2 ). This is particularly important for heavy quark physics since, in 
current lattice simulations, the typical heavy quark mass tuq is rather large, niQa ~ 0.3-0.6. 

Since the coefficient of the Clover term is known non-perturbatively, the hadron spectrum 
is definitively improved to 0(a 2 ). Unfortunately, the program of removing all the 0(a) 
corrections in the operator matrix elements out of the chiral limit has not been completed 
yet, although strategies to this purpose already exist || [R]]. For this reason, in some cases, 
we have used the improvement coefficients (1>a > cy, b m ) evaluated at first order in (boosted) 



perturbation theory [11], thus leaving us with 0{p? s am) corrections, where m is the relevant 
quark mass. 

After a careful analysis of the systematic uncertainties present in the extraction of the 
physical results, by assuming quite conservative errors, and bearing in mind the systematic 
effects due to the quenched approximation, the main results of our investigation are the 
following: 

i) For D mesons we find: 

f D = 211 ± 14 1° 2 MeV, f Ds =231 ± 12 ±g MeV, ( = 1.10(2), 



D 



fo- = 245 ± 20 lg MeV, f D * = 272 ± 16 1° MeV, =1.11(3), (1) 

where /d* and fr>* are the vector-meson decay constants. The latter quantities are not 
measured experimentally, but enter the calculation of two-body non-leptonic -B-decays 
computed using factorization [12|. Thus, they are useful for checking the factorization 
hypothesis with charmed vector mesons in the final states. 

ii) For B mesons, we find: 

f B = 179 ±18 tf MeV, f Bs = 204 ± 16 Ijf MeV, 0^=1.14(3)1?, 
f B , = 196 ± 24 l^ 1 MeV, f B * = 229 ± 20 1^ MeV, = 1.17(4) 1°. (2) 
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Following Ref. 



1 13], we have also directly computed the ratio 



■y^- =0.78 ±0.04 tj 1 (3) 

IDs 

from which, using the experimental result, f^ p ) = 254 ± 31 MeV QQ, we find 

fB = j^x = (198 ± 24(exp.) ±f (theo.)) MeV , (4) 

JD a 

Although with a larger error, the result in Eq. (||) is well compatible with the value 
given in (g). 

iii) To reduce the effects of the quenched approximation, we have also used the Grinstein 
double-ratio method (illustrated below), obtaining 



td = ^ = 1-19(5) , r B = ^ = 1.23(6) , ^- = 0.71 ± 0.04+*° . (5) 

JD IB IDs 

The latter ratio would give as the best estimate for fs- 

f B = (180 ± 26(exp.) j3jj(theo.)) MeV . (6) 
With the double ratio method we also obtained 

r B /r D = 1.03(4) . (7) 

iv) We made a detailed study of the hyperfine splitting and of the scaling laws for masses 
and decay constants, as predicted by the heavy quark symmetry. The results of this 
study can be found below. 

We now give the details of our analysis and of the methods used to extract the different 
physical quantities. Since most of the techniques are by now standard and have been de- 
scribed ad abundantiam in the literature^, we only focus on those points which are either less 
common or new. More details on the calibration of the lattice spacing and on the extraction 



of the hadron masses and matrix elements can be found in Refs. [17, jig] . 

The remainder of this paper is organized as follows. In Sec. ||, we list the main parameters 
of our simulation and introduce the basic notation necessary for the discussion of the results. 
The heavy-light meson masses and decay constants in lattice units are also given in this 
section. Since the systematic effects related to the extrapolation/interpolation to the physical 
point, although related, are quite different in the two cases, we present separately the physical 
predictions for D and B mesons, in Sec. ^ and Sec. ||, respectively. In Sec. ||, we discuss 
the scaling laws predicted by the heavy quark effective theory (HQET) and other related 
subjects. 



1 This value has been recently updated by the same authors and reported to us by F. Parodi. 
2 Reviews, with complete lists of references, can be found in Hq]. 
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2 Lattice results 



In this section, we give the essential information about our numerical calculation and establish 
the basic notation. We then present our results for the heavy-light meson masses and decay 
constants in lattice units. 

The numerical simulation has been performed on a 24 3 x 64 lattice, at = 6.2, in the 
quenched approximation. All results and errors have been obtained with a statistical sample 
of 100 independent gauge field configurations, using the jackknife method with different 
decimations. We have used the non-perturbatively improved lattice Clover action, with 
c sw = 1.614 ||]. We work with four values of nu g ht, and four Kheavy'- 

• 0.1352 (k £i ); 0.1349 (« /a ); 0.1344 0.1333 (k u ), 

• 0.1250 (K hl ); 0.1220 0.1190 (^ 3 ); 0.1150 ( Kh4 ). 
From the study of the light-hadron spectrum, we obtained 

- a-^mK*) = 2.75(17) GeV, 

- K cnt = 0.135845(25), 

- K q = 0.135804(26) , 

- k s = 0.13482(12) , 

where K q , corresponds to the light quark mass m q (with q = u,d), and k s to the strange- 
quark mass, m s . The above values have been obtained from the physical pion and kaon 



masses, by using the method of physical lattice planes [19|. All details regarding light hadron 
spectroscopy and decay constants, can be found in Refs. |l^, |lS[ |. 

For the mass spectrum, following the standard procedure, we measured suitable two-point 
correlation functions, from which we can isolate the lowest lying states 



C,At) = £<0|J(*, i)jt(0)|0> m -^e- MrT/2 ™sh [ Mj (I - t) 



(8) 



where J = Jps = Q^fsq, or J = Jy = Q"f k q. In Fig. |], we show the effective masses for 
the pseudoscalar and vector heavy-light mesons at fixed heavy quark mass. By inspection, 
we established the fit intervals t £ [20,28], and t £ [22,28], for the pseudoscalar and vector 
cases, respectively. The resulting pseudoscalar and vector masses in lattice units, as well as 
the matrix elements, Z PS = \ (PS(p = 0)|J PS |0)| 2 and Z v = \(V(p = 0; A)| Jy\0}\ 2 , are listed 
in Tab. § 
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Figure 1: Effective masses of heavy-light pseudoscalar and vector mesons as a function of the time 
in lattice units. In each figure, the heavy quark mass (corresponding to = 0.1220,) is fixed, and 
combined with four different light quark masses. 
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Table 1: Mass spectrum of heavy-light pseudoscalar and vector mesons in lattice units. 

We used the standard procedure to extract the pseudoscalar and vector decay constants. 
This procedure consists in calculating the ratios 



E s (p&t)pm ~ Wzp v v 2 

^(VifrWm ^ M*F*e- M v% cosh (V y (|-t)) (10) 

where we assumed the usual definitions 

(Q\A o \PS{p = 0)) = iFpsMps, and 

(0\V i \V(p = 0;X)} = ie\ X) F v M v . (11) 

We denote decay constants and meson masses in lattice units by capital letters, and the hat 
reminds us that the quantity is improved and renormalized. In practice, one first partially 
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improves the bare lattice currents (for clarity, we write the lattice spacing a explicitly): 



(0\A \PS(p = 0)> -» (0\A \PS) + c A (0\ad P\PS) = iM PS (Fj® + c A aF^ s ), 

(0\V l \V(p = 0;X)) - (0|^|y) + cy(0|a9or j0 |y) =iM y ef ) (F^ 0) +c y a4 1) ), (12) 
and then multiplies the currents by suitable overall factors 



F P = Za (1 + ^am) (F> Uj + c A aF ( p>), [Z A (m) =Z A (1 + b A am)] 

4° > + cyaF«), 



Fy = Zy (1 + oyam) (F^ + cyaF^), [Zy(m) = Z A (1 + b v am)] . (13) 



In the calculation of the different correlations above, when the lowest state is well isolated, 
we may use 



aF^ s = ^^- S mh(Mps), and aF« = sinh(My). (15) 

The values of the decay constants are given in Tab. 

The improvement coefficients and the renormalization constants are catalogued in Tab. 3, 
where we also display the one-loop results obtained by using boosted perturbation theory 
(BPT) at /3 = 6.2 [11, ^O|0. Recall that the corrective coefficients bj enter with the "average" 



quark mass defined as am = arriij = \ (ami + amj), where the bare mass is the one derived 
from the vector Ward identity, namely 

amt = -(—- —-) . (16) 



2 y K% ^crit 

In the following, we denote by m q and toq the generic light and heavy quark masses, whereas 
the quark masses expressed in terms of the corresponding hopping parameters, as in Eq. (|l~6l), 
are denoted by mg or m^. 

Note that, in spite of the non-perturbative determination of cy, we used the perturba- 
tive value Cy PT = —0.026. Firstly, we find the non-perturbative result of Ref. H, Cy P = 
—0.214(74), surprising because it is one order of magnitude larger than c v PT . This possibil- 
ity is not excluded a priori, but it is difficult to accommodate it in the pattern of all other 
improvement coefficients: when known non-perturbatively, their value is always close to the 
corresponding (boosted) perturbative one and never differs by one order of magnitude. Sec- 
ondly, by using c v p = —0.214(74), the ratio of the vector to the pseudoscalar meson decay 
constants /#*/ ' fjj badly fails in approaching one, as Mh increases, contrary to what is pre- 
dicted by heavy quark symmetry. More details on this scaling law will be given in Sec. H For 
these reasons, we find it safer to use the c v PT . We believe that the preliminary determination 
of c v p in Ref. || has some problem and prefer to wait for the final results. 



3 This corresponds to the use of g 2 — 1.256, in the perturbative formulae. 
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Table 2: Heavy-light decay constants in lattice units. 



3 D-meson spectrum and decay constants 

In this section, we discuss the D-meson spectrum and decay constants. Preliminary results 
of this study were given in Ref. |^][]. 

In Tab.|4], we tabulate the results for the heavy-light meson masses, M^(m/j, rri£ ), obtained 
from a linear extrapolation (interpolation) in the light quark mass (to reach q = u,d and s). 
This was achieved by using the method of physical lattice planes. In Ref. fi~7|| , we extracted 
the M n , and the hypothetical pseudoscalar M Vss , which (when squared) are proportional 
to m q and m s , respectively. For the generic physical quantity in the heavy- light meson 
sector $if(mh,me), we use the following form of fit 

$H(m h ,m ei ) = a h + f3 h M PS (m et ,m £i ) + j h {M 2 p S {m lv m^)) 2 , (17) 
4 See [^D for preliminary results from the UKQCD collaboration and the APETOV group. 
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Renormalization constants 
(in the chiral limit) 


Quantity 
BPT 
Non — perturbative 


Z v 
0.846 
0.793 


Z A 
0.862 
0.809 


Coefficients for the Improvement 
of the Bare Operators 


Quantity 
BPT 
Non-perturbative 


Cy 
-0.026 

-0.214(74) 


CA 
-0.012 
-0.037 


Coefficients for the Renormalization Constants 
Improvement (due to the Explicit Mass Term) 


Quantity 
BPT 
Non-perturbative 


bv 
1.242 
1.404 


b A 
1.240 

NOT CALC. 



Table 3: Improvement coefficients. In boosted perturbation theory g 2 = 1.256. For the perturbative 
Zj's, we used c sw = 1.614. The values which have been used in our numerical calculations are 
marked in bold. 



where the heavy quark mass (i.e. K/J is kept fixed. The coefficients of such a fit, ah, 0h, 
jh, are then used to obtain 5f/(m/,,m 5 ) and Sjj^m/i, m s ), by inserting on the r.h.s. of fllT]), 
M 2 and M 2 , respectively. In practice, it turns out that the linear (7^ = 0) and quadratic 
(lh 0) Hts give essentially the same results for any physical quantity considered in this 
study []. In Fig. ^, we show this effect for the pseudoscalar decay constant. Therefore, in all 
results that we present in what follows, whenever a quantity with light quark flavor q and/or 
s is mentioned, it means that the linear fit in Eq. (|i~7| ) is performed, i.e. = 0. 

Having fixed the light quark mass, we now want to interpolate in the heavy quark mass. 
In the framework of the HQET, the functional dependence of Mjj on ttjq is: 

M «=" ! «{ 1+ 4 + ^5 (Al+ * A2)+0( 4 ) } (18) 

where A is the so-called binding energy, Ai_2 are the first (flavor-spin) symmetry break- 
ing corrections (describing the kinetic and chromomagnetic energy), and k = 3(— 1), for 
J p = 0~(1~). The improvement of the quark mass brings in the quadratic terms in m/j, 
i.e. mh — ► rrihil + b m mh), and distort all the coefficients in the expansion (|l8|). The term 
of order m? originates only from the lattice artifacts, and thus is always proportional to 
b m . The interplay between power corrections in 1/mg and discretization effects, however, 

5 As expected, the results obtained from a quadratic fit inflates the errors in extrapolated results. 
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0.9721(61) 


0.0870(36) 
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0.6409(35) 
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0.0896(51) 



Table 4: Mass spectrum and decay constants of heavy-light pseudoscalar and vector mesons. hi~s 
and hi-q denote mesons composed by a heavy quark with mass and a strange or a light (u, d) 
quark, respectively. All the results are expressed in lattice units. 



modifies the "effective" value of b m , i.e. the coefficient of the quadratic term in m^. To 
investigate this point, we study the behaviour of Mh in mh, at fixed light quark mass, m q . 
In the D case, we use the following expression: 

M H (m hi ,m q ) -am D = A(m hi - m charm ) [1 + B (m hi +m charm )) (19) 

where amp is the experimental meson mass in lattice units, amp = 0.68(4) (similarly we 
fit Mh* (m^, m q) ~ amn*, etc.). From the fit of our data to Eq. (|l9|), it turns out that the 
resulting value for n c harm is stable for B £ [—0.4,-0.2]. The minimum x 2 is reached for 
B = —0.32. We have also performed the linear fit (corresponding to B = 0), and the fit with 
B = b m = -0.652 0. The different values that we obtain for K c harm with different fits 
(linear, quadratic or using b m ) differ by about one per mille. We quote 

Kcharm = 0.1231(14). (20) 

It can be argued that a fit of the spin-average mass Mh = (3My + Mps)/4, to extract 
^charm is more suitable, because spin forces of 0(1/ mQ) are canceled in this combination 
(see (|l8|)). For B = —0.32, corresponding also in this case to the minimum x 2 , we obtain 
Kcharm = 0.1232(14). Since the differences for the D— meson masses and decay constants as 
obtained by using the two values of n c harm is very small, in the following, whenever we refer 
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Figure 2: Fits of Fp$ in the light quark mass, at fixed rrih- The filled circles denote data directly 
measured. The dashed curve and empty circles refer to the linear fit and extrapolated points. The 
dotted curve and empty squares to the quadratic fit and extrapolated points. The heavy quark mass 
corresponds to Kh 2 - 



to Kcharmi the value (|^) is understood. Using Mr, as a physical input (to fix K c h arm ), we 
can make several predictions for other meson masses 



M D = input ; M D * = 0.725(42), and 
M Ds = 0.733(46) ; M D * = 0.768(45). 

which in physical units give 

m D * = 1.992(24) GeV , m Ds = 2.013(18) GeV , m D * = 2.110(21) GeV 

to be compared to the experimental numbers 



rn 



(exp.) 
D* 



2.008 GeV; mgf 



1.968 GeV; 



m 



(exp.) 



2.112(27) GeV . 



We obviously fail to obtain the experimentally measured mass-difference. We get 

m D *-m Ds = (97 ± 12) MeV. 
which is to be compared to (mp* — mo s )^ exp ^ = 143.8(4) MeV. 
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An alternative procedure is to consider the ratio C vv /C pp , from which the vector- 
pseudoscalar mass difference can be directly extracted. By using this method we get 

M D * -M d = 0.0354(39) -» m D * - m D = 97(15) MeV, 

M D * - M Ds = 0.0354(30) -» m D% - m Ds = 97(13) MeV; 

which confirms that the ('spin') mass difference is systematically smaller than the experi- 
mental one, regardless of the procedure we use. Since we found a reasonable agreement for 



the hyperfine splitting in the light-quark sector [17], the discrepancy in the heavy-quark case 
is probably a signal of large residual 0(a 2 ) errors. We believe that the discrepancy cannot 
be entirely attributed to the use of the quenched approximation ^. 

Our results for the hyperfine splitting are shown in Fig. ||. From that figure, we note 
the qualitative agreement between the dependence of the splitting on the light-quark mass 
measured on the lattice and its experimental counterpart. Moreover, the dependence of the 
hyperfine splitting on the meson mass is not dramatically larger than the experimental one, 
represented by a gray line in the figure. This is to be contrasted to the case of the unimproved 
Wilson action, where the lattice slope is by far larger than in the present case [pi]], showing 
a clear effect of improvement, although insufficient to describe the experimental numbers. In 
Tab. |, we list the results extrapolated in the light quark mass, at fixed m^. 



0.08 



K h =0.125 

" 3 0.06 - K h =0.122 



k. =0.1 19 

*J*0.04 - 1 



0.02 



0.6 0.7 0.8 0.9 1.0 1.1 

(3M Vm+ M PSu )/4 

Figure 3: Hyperfine splitting for heavy-light mesons in lattice units. The gray line shows the 
( approximate ) experimental slope. The stars mark the physical m 2 D , —m 2 D and m? D , —m 2 D splittings. 



'Quenching is always a joker-argument when we are unable to solve or explain a problem in lattice calculations. 
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heavy 


q 


s 


hi 


0.393(53) 


0.406(50) 


h 2 


0.367(57) 


0.383(52) 


h 3 


0.342(61) 


0.361(54) 


/14 


0.310(68) 


0.330(57) 



Table 5: Hyperfine splitting, Mjj t — Mfj in GeV 2 



We now discuss the physical results for the D-meson decay constants. We first extrapolate 
the results, in the light quark mass, at fixed m^, by using a formula similar to Eq. (|17|), i.e. 



F H (m h ,m e J = a F h + /?£ M PS (m^ , m u ) . 



(25) 



The results of this extrapolation, m,£. — ► m q and — > m s , are given in Tab. |3]. 

To handle the problem of extrapolation in the heavy-quark mass, at fixed light-quark 
mass, we rely on the heavy quark symmetry. The relevant scaling law is 



fH 



$'(m g ) 
§{m H )m H 



1 ! ...j , (26) 

where $(m#), 3>'(m#) depend logarithmically on the mass, e.g. $(m#) ~ a s 2 / 6o (m#)(l + 
C(a s )) []. In the interval of masses considered in this study, the logarithmic corrections are 
negligible. For this reason, in our fits, we used 



$1 



$ 2 



(27) 



where $o; $1 an d ^2 are constants. At the physical point Mjj = amp (corresponding to 
Kh = K-charm), we read off the value Fp> in lattice units. To express it in physical units, 
one simply multiplies by a~ 1 . The same procedure can be used for the vector-meson decay 
constants. 

Another possibility, is to consider the ratios RH(mh,rri£) = Fjj(mh,me) / 'Fps(m q ,me) 
and Rh* = Fij*(mh,mt)/Fy(rn q ,rn£), and to extrapolate Ruim^mi) in mi and by 
using eqs. ( p5| ) and (^7j), with the obvious replacement Fjj — > Rh (F ntP = Fpsy(m q ,m q ), 
Fr,k* = Fpsy(m q ,m s )). The physical values of the decay constants are then obtained by 
using 



Sd = RH(m charm ,m q ) x 



(exp.) 



fD 3 



RH a {m charm ,m s ) x f^ xp,) 



(28) 



7 We prefer to give the scaling law in terms of the hadron mass m# rather than the heavy quark mass. 
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and similarly for the vector mesons 

Id* = Rh* (m charm , m q ) x ff xv ^ , f D * = R H * {m charm , m s ) x f^T ] . (29) 

The experimental values of the decay constants that we use are the following ones f24|| : 
j(exp.) = 131MeVj fi^P-) _ i 6 o MeV, f { p Xp - ] = 208(2) MeV, f^7' ] = 214(7) MeV. 

The results are given in Tab. |6|. We also give the decay constants obtained by including 
the KLM factor which we discuss in Sec. ||. The differences can be used for an estimate of the 
residual 0(a 2 ), discretization errors in the determination of the matrix elements. In Tab. 0, 
we list the corresponding results for vector mesons. 



No KLM factor 


fn q [MeV] 


f Ds [MeV] 


linear in 1/M# quad, in 


linear in 1/Mjj g quad, in 1/Mh s 


fot using F n 
f De using F K 
f De using aT l (m K *) 


201(22) 200(21) 
213(14) 212(15) 


239(18) 238(16) 
233(11) 232(12) 


With KLM factor 


fo q [MeV] 


Id. [MeV] 


linear in 1/Mh quad, in 1/Mh 


linear in 1/Mh s quad, in 1/Mh s 


f De using F n 
f D( using F K 
fo t using a~ l (m K *) 


199(22) 198(21) 
211(14) 210(15) 


237(17) 236(16) 
231(12) 230(13) 



Table 6: Pseudoscalar decay constants for D-mesons using the scaling law in Eq. (|27D. Results 
including the KLM factor discussed in the text, are given in the lower part of the table. 



Whether we use a linear or a quadratic fit to interpolate to n c harm, our results in the 
.D-sector remain practically unchanged. In order to illustrate the stability of the results for 
D mesons, we also show in Fig. |||, the results of the linear and quadratic fits in 1/Mh- 

The observed stability makes these results quite remarkable: we use the non-perturbatively 
improved action; the operators and the renormalization constants are also improved; the re- 
sults obtained by using the heavy quark scaling laws are unchanged, regardless of whether we 
take quadratic (1/M H ) corrections into account or not; the results are practically insensitive 
to the presence of KLM factors; there is no important dependence on the quantity chosen 
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,M „ 0.09 - T 

M . I • H mesons 



O H s mesons 
— Linear fit 

0.08 h Quadratic fit 



0.07 



0.06 - 



i I 



7.00 1.20 1.40 1.60 

1/M H 

H (q,s) 



S 0.11 



"„ 0.70 



0.09 - 



0.07 



0.07 



7.00 



7.20 



• H* mesons 
O H* s mesons 
— Linear fit 
Quadratic fit 




1.40 



1/M H , 



7.60 



(q,s) 



Figure 4: Results of the linear and quadratic fits for pseudoscalar and vector mesons, in lattice 
units. The gray dashed lines correspond to fu (foj, and fo* (Id*) respectively. 
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fD* 


[MeV] 


Id; 


[MeV] 


No KLM factor 


linear in 1/Mh 


quad, in 1/Mh 


linear in 1/Mh 3 


quad, in 1/Mh s 


f D * using F p 


246(30) 


244(32) 


- 


— 


f D * using F K * 


— 


— 


255(17) 


253(18) 


f D * using aT l (m K *) 


248(19) 


246(21) 


275(15) 


273(16) 




fD* 


[MeV] 


/d* 


[MeV] 


With KLM factor 


linear in 1/Mh 


quad, in 1/Mh 


linear in 1/Mh 3 


quad, in 1/Mjj s 


Id* using F p 


243(31) 


241(33) 






Id* using F K * 






252(18) 


251(19) 


f D * using a~ l {m K *) 


245(20) 


243(21) 


272(16) 


270(16) 



Table 7: Vector decay constants for charmed-mesons obtained by using the scaling law in Eq. (p7j). 
Results, including the KLM factor discussed in the text, are also given. 



to fix the physical normalization. The errors that we quoted in (Q), are obtained in the 
following way: a central value is fixed by the result obtained from the linear fit in 1/Mh, 
with the scale fixed by mjf«, and the KLM factor included; we quote the statistical error 
as estimated using the jackknife procedure; all the residual differences are lumped into the 
systematic uncertainty (the difference between the central values of the results obtained by 
using different quantities for the scale fixing, and the difference with the central value of the 
result obtained from the quadratic fit in 1/Mh)- It is also worth noticing the remarkable 
stability of the ratio /d s //z> (see Eq. @), although it may be questioned whether we are 
really able to predict the SU(3) breaking properly in the quenched approximation. More 
discussion on this point will be given in Sec. ||. 

4 B mesons 

In this section, we present the results of the extrapolation of the decay constants to the B 
mesons, and discuss the discretization errors in the extrapolation. 

When extrapolating the raw data obtained for m/j ~ m c h arrn to the -B-sector, two im- 
portant effects may arise. On the one hand, the inclusion of the quadratic term ^2/Mfj in 
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Eq. (p7D may change appreciably the results of the linear fit, on the other the 0{a) correc- 
tions (0(a) terms proportional to cy,A and by, a in (|l3|)) become much larger. This is to be 
contrasted to the case of -D-mesons, where the inclusion of the quadratic corrections leaves 
the results essentially unchanged, cf. Tabs. |6| and [?]. 

The effect of ca, cy, by and 6 a is sizable for the scaling behaviour of f B)B *. Note also 
that if we used c v p , this effect would be huge for the vector decay constant. For instance, 
in the range of quark masses considered in our simulation, the renormalization constants 
Zy^irn), defined in Eq. (|13"D, increase by 20 4- 50%, relatively to their values in the chiral 
limit. Since Zy(m) and Za{jti) are multiplicative factors, their effect is very important for 
the extrapolation to the 5-sector. This is illustrated in Fig. H when Za{jtl) is included, we 
note that the quadratic fit is more desired, although the linear one is compatible with the 
data. The embarrassing point is that the values of f B and f B * , as obtained from the linear 
and quadratic fit, are hardly compatible, see Tabs. || and ||. This is particularly pronounced 
for Bs mesons. 





fB q 


[MeV] 




[MeV] 


No KLM factor 


linear in 1/Mh 


quad, in 1/Mjj 


linear in 1/Mh s 


quad, in \/Mh b 


f Be using 


176(19) 


208(27) 






f Bl using F K 






217(14) 


255(20) 


f Be using a~ l {m K *) 


187(19) 


220(25) 


212(16) 


249(20) 




fB q 


[MeV] 


Ib s 


[MeV] 


With KLM factor 


linear in 1/Mjj 


quad, in 1/Mh 


linear in 1/Mh s 


quad, in l/M Bs 


f Bg using F n 


170(18) 


193(25) 






f Be using F K 






209(13) 


238(19) 


f Be using a' l (m K *) 


179(18) 


205(24) 


204(16) 


232(19) 



Table 8: Pseudoscalar decay constants for B-mesons using the scaling law in Eq. (P?|). Results 
with the KLM factor included are listed in lower part of the table. 



The curvature in the fit to f B could partially be induced by 0(a 2 ) terms, still present in 
the calculation of the matrix elements. A possible way to account for some of these effects 
is through the so-called KLM factor [pfif| . In our case, this means that, besides the factor 
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0.10 



0.09 - • F h M h' /2 ( bare ) 



0.08 - 

0.07 - 
0.06 

0.05 - 



0.10 



0.08 - 

0.07 - 

0.06 - 

0.05 - 



I 



q q 



\ i i * 



1.00 1.20 1.40 1.60 



1/M H 



1/2 



0.09 - OF n q M H q ( no KLM) 

• F H M H " 2 (with KLM) 



q q 



7.00 7.20 7.40 7.60 

1/M 



Figure 5: Heavy-light pseudoscalar decay constant as a function of 1/Mh- The two figures illus- 
trate the influence of the renormalization constant ZAijn) on the 1/Mh dependence of the decay 
constant. 
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No KLM factor 


f B * [MeV] 


fB* [MeV] 


linear in 1/M B quad, in 1/M B 


linear in l/M Bs quad, in l/M Bs 


f B * using F p 
f B * using F K * 
f B * using a~ x (mK*) 


204(34) 239(39) 
205(25) 241(32) 


222(22) 260(25) 
239(21) 280(24) 


With KLM factor 


f B * [MeV] 


f B * [MeV] 


linear in 1/M B quad, in 1/M B 


linear in l/M Ba quad, in l/M Ba 


f B * using F p 
fs* using F K * 
f B * using a^irriK*) 


194(32) 225(37) 
196(24) 227(30) 


213(22) 241(24) 
229(20) 260(23) 



Table 9: Vector decay constants for B-mesons using the scaling law in Eq. fl2T|). Results including 
the KLM factor discussed in the text, are also given. 



(1 + bjma) already included in the definition of the renormalized currents fll3|), we may try 
to include the effects of higher order terms in ma, by using the following relation 



Zj(m h ,m e ) 



Zj(0) 



y/1 + amh y/1 + ami 
1 + am 



(1 + bj am) 



~ Zj(0) (1+ bjam) + 0(a 2 



(30) 



where am = 
Equation ((3 1 



(amh + o,m^)/2, and ami is the usual expression for the bare quark mass dig 

\/T+ am g 



is a consequence of the redefinition of a quark field, q — > V i + am g (m 
the KLM way), which comes from the comparison of the /ree lattice quark propagator to 
its continuum counterpart. The results which include the KLM correction are given in the 
lower part of Tabs. [], [?], || and |9[ In the case of .D-mesons, the effect of KLM is indeed 
negligible. In the case of B-mesons, we observe a slight change in the central values, e.g. 
f B = 187 MeV — ► 179 MeV. However, the distance between the values obtained with linear 
and with quadratic fits remains essentially unchanged. In the absence of a larger range of 
masses, we are unable to reduce the difference between results obtained with the linear and 
quadratic fits. As it has been done for D-mesons, we quote the results of the linear fits as our 
central values, and include in the systematic error the differences between our central values 
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and i) the results from the quadratic fit; ii) the results without the KLM factor incorporated; 
in) the results obtained by using other quantities (fiCi fn) to extract the physical values. Our 
final results are those given in Eq. (||). 

5 Scaling laws and related issues 

In this section, we discuss several interesting quantities for the study of the scaling laws 
predicted by the HQET, and their validity in the range of quark masses between the charmed 
and the bottom one. We introduce several ratios of decay constants which are useful to get 
some physical information. 

We first consider the scaling law for the decay constants. The results for the coefficients 
in Eq. (|27|), as obtained from our fits, are given in Tab. 10. To translate these coefficients 
into physical units, we have used a~ 1 (mx*)- The leading term from the linear fit, $o = 



fit 

parameters 


mi = m q 
linear quadratic 


mi = m s 
linear quadratic 


$£ s [GeV 3 / 2 ] 
$f s /^ s [GeV] 
^/W [GeV] 


0.48(5) 0.66(11) 
-0.75(6) -1.60(22) 
1.03(8) 


0.56(5) 0.74(8) 
-0.83(5) -1.70(16) 
1.08(6) 


<S>% [GeV 3 / 2 ] 
<Dj7*o y [GeV] 
yj^l^ [GeV] 


0.51(7) 0.70(12) 
-0.63(9) -1.62(24) 
1.09(9) 


0.61(6) 0.81(10) 
-0.74(6) -1.65(17) 
1.08(6) 



Table 10: Fit parameters in physical units for pseudoscalar (PS) and vector (V) heavy-light 
mesons. 



0.48(5) GeV 3 ' 2 , is in good agreement with the findings of previous studies [13, 26, 27, 28] 
We also note that this value is compatible with the results of QCD sum rules 
(0.4 -7- 0.6) GeV 3 / 2 , when the large perturbative QCD corrections are included 0. 



9 , $r 



Following the Refs. [21 



we now consider the spin scaling relation on the lattice: 



U(M H ) J " 



fH* 



i0+ w~ H + w H 



(31) 



Without these corrections, the result would be <&o = 0.30(5) [E9| 
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where Mh = (3M#* + M#)/4 is the spin averaged mass (which we already used in Sec. |3|), 
and £0,1,2 are parameters which we obtain by fitting our data. From heavy quark symmetry, 
one expects that £0 = 1 (up to logarithmic corrections). For completeness, we tabulated M h 
and ftf/fx* in Tab. |ll|. The results of our fits in physical units, are 



'heavy flavor' 


K hl K h2 K h3 K h4 


M H = (3M H * +M H )/A 
M Hs = (3M H *+M Hs )/4 


0.6705(48) 0.7620(53) 0.8486(56) 0.9579(61) 
0.7028(58) 0.7940(58) 0.8802(58) 0.9890(60) 


U(M H ) = fulfil* 
U(M Hs ) = f Ha /f H1 


0.851(30) 0.869(29) 0.879(29) 0.888(31) 
0.844(22) 0.858(21) 0.867(21) 0.876(22) 



Table 11: Spin averaged masses and ratios of pseudoscalar and vector decay constants. For Mh, 
the light quark mass q = u,d, is understood. 



(lin.) Co = 0.997(68) , £i/£o = -0.23(11) GeV , 

{quad.) Co = 0.89(12), £i/£o = 0.17(49) GeV , x/&/£o~ = -0.67(18) GeV , (32) 

where the physical values were obtained by using a _1 (m^»). Data points, and extrapolated 
values, are displayed in Fig. ^. We see that the scaling law is very well satisfied by using the 
linear fit. The inclusion of the quadratic term, although irrelevant in the directly accessible 
region of the meson masses, produces large deviation from the expected extrapolated value 
£0 = 1, as Mh — > 00. Thus, by using the linear fit, we arrive at 

U(M D ) = 0.860(28), C/(M B ) = 0.933(47), (33) 

and 

U(M Ds ) = 0.868(21), U(M Bs ) = 0.915(33). (34) 

We end this section by presenting a set of ratios which may be explored in order to extract 
the physical decay constants by using some measured quantities. 

• As it was suggested in Ref. the decay constants can be conveniently presented in 
terms of fo s , which is already measured 9 : 

= 0.78 ±0.04 tl 1 ] ^ = 0.88 ± 0.03 1 2 ; 

JD S ' IDs 
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_F 1.1 - 



1/M„ 



1/M„ 



7.0 - 



OS 



OA 



0' 



0.7 



0.0 



0.5 



7.0 



7.5 
1/M 



H 



Figure 6: Ratio of the heavy-light decay constants. The linear fit approaches very well the expected 
asymptotic value U(Mh — > oo) = £ = 1. The results refer to mesons with the light quark 
extrapolated to q = u, d. 



^ = 1.06 ±0.05; ^p- 

JDs ID a 



1.17 ±0.03; 



^ = 0.85 ±0.07 ^ = 0.99 ±0.05 1^ 4 . 

JD S IDs 



(35) 



The error estimates are obtained in the same way as in Sec. |j. /b//d s is the value 
which has been used in Eq. (|j), as an alternative way to extract the value for fs- 

Other phenomenologically interesting ratios for testing the factorization hypothesis in 
non-leptonic modes, are (see [pCj] ): 



1.29(14); 



fDs 



1.66(19). 



In the quenched approximation, the S£/(3)-breaking parameter rx — 1 = fx/fn 
is expected to be smaller than its experimental value. A smaller value of tk 



(36) 

- 1, 
1 

is predicted by one- loop quenched chiral perturbation theory [31|, and is verified in 
numerous simulations (with either unimproved or improved actions and operators |HT 
A similar effect is also expected for r H = fujfn (r D = /d s //d, r B = /b s //b) 
In this respect (in the hope of reducing the quenching errors), it may be interesting to 
examine the Grinstein-type double ratio 9\h = rff/ro p3|- From our data, we have 



m Hhi = 0.995(3) , K Hh2 = 1.003(4) , ^ = 1.009(6) , 9^ = 1.014(9) , (37) 
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which upon an extrapolation to the B-meson mass, amounts to 



m (Un) = 1-0 35(i7) ; y^ uad ) = 1.028(33) . (38) 

Using and r D = 1.10(2) from (|), we have vb = 1.134(34), in perfect agreement 

with the direct determination, given in (g). 

The double ratio can be used to estimate the quenching errors in the predicted values 
of T£> and rs- To this purpose, we define 



(exp.) 

ru -(f) . (39) 



r H ( Ik\ 
r K 



where r# and rx, are obtained in the quenched lattice calculation. Using our data 
(r K = 1.12(5) g|) and (f K /U)( exp ' ) = 1-22 H, we end up with 



^HnMguad.) = hlg ^. (4Q) 



-(lin.&tquad.) 



1.23(6), (41) 



which are ~ 9% larger than the results obtained directly and quoted in Eqs. (|T]) and (|2|). 
If this difference of 9% is the realistic estimate of the quenching errors, they are much 
smaller than the pessimistic estimate of Ref. |3^], where ~ 20% of (quenching) error 
was predicted. Note that the ratios fjj and fs, do not depend on the fit we use (linear 
or quadratic) . A similar game with Jb / Jd 3 results in 

■p- = 0.71 ±0.04 +1° (42) 

JDs 

which gives 

Sb = (180 ± 26(exp.) ^(theo.)) MeV, where we accounted for the experi- 
mental value for f^ xp '\ This result agrees with the value we reported in (||). 
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6 Conclusion 



We have analyzed masses and decay constants of heavy-light pseudoscalar and vector mesons, 
using the non-perturbatively improved action and currents. Particular attention has been 
paid to the errors coming from the extrapolation in the light and heavy quark masses. 

We find that the hyperfine splitting is definitely below the experimental value, in spite of 
the improved action. 

The values predicted for the decay constants of D mesons are extremely stable against 
variations of the fitting procedure, inclusion of KLM factors etc. Thus we believe that the 
main error on these quantities is the quenching error. 

On the contrary, we find larger uncertainties for the S-meson decay constants, mainly 
due to the amplification of discretization effects when extrapolating to the 6-quark mass, 
and to the uncertainty in the extrapolation procedure. In spite of these uncertainties, and 
of the fact that our results are obtained at a single value of the lattice spacing, we believe a 
value of Jb much lower than 170 MeV rather unlikely. Indeed, for (3 > 6.0, with Wilson-like 
fermions at fixed lattice spacing, almost all lattice calculations give fs larger than 160 MeV. 
This value has been quoted as the "world average" obtained in Ref. [34], after combining 
results obtained with propagating quarks, with those obtained using some effective theory, 
as NRQCD |$5|, or the FNAL action [36]. Low values of fs with propagating quarks are 
obtained only after extrapolating in a to the continuum limit [^, [?4|], with procedures which 
we believe are questionable (for example by including data at low values of j3, i.e. too 
close to the strong coupling regime). Our results, which should have smaller discretization 
errors than other calculations at fixed lattice spacing, confirm a value of fs (in the quenched 
approximation) larger than 170 MeV. A (rather) indirect evidence that a larger value of fs 
is preferred can be obtained by combining Jb/ Jd 3 from the lattice with the experimental 
value of fo s - This gives fs — 180 -j- 190 MeV, with an error of about 40 MeV. Finally, we 
used the Grinstein double-ratio method, in order to try to reduce the quenching errors for 
(ratios of) decay constants. 
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